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ABSTRACT 


The  compression  response  of  thick  (6.35  mm)  and 
thin  (1.60  mm)  (0./90)  AS4/3501-6  laminates  was 

compared  for  material  in  the  dry  and  moisture 
saturated  conditions.  The  thick  specimens  failed  at 
1055  MPa  dry  and  740  MPa  wet.  Comparable  strengths 
for  the  thin  specimens  were  1090  MPa  and  957  MPa. 

There  was  no  statistical  difference  in  strength  for 
the  dry  material  of  either  thickness.  The  thick 
specimens  absorbed  significantly  more  water  than  the 
thin  specimens,  and  consequently,  showed  a  greater 
loss  in  compression  strength. 

The  elastic  modulus  was  not  affected  by  moisture 
in  the  thick  or  thin  specimens. 

The  gross  failure  mode  of  the  dry  and  wet 
specimens  of  either  thickness  was  very  similar  with 
slightly  more  delamination  occurring  in  the  wet 
specimens.  The  thick  specimens  typically  exhibited 
angular  cracks  of  27°  to  30°  partially  across  the 
specimen  thickness.  This  was  not  observed  in  the  thin 
specimens. 

Finite  element  techniques  were  used  to  determine 
shear  stresses  in  the  matrix  for  a  fiber-matrix 
composite  plate  with  an  initial  imperfection,  and  loaded 
in  compression.  These  matrix  shear  stresses  were  used 
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in  a  compression  failure  theory  which  showed  reasonable 
agreement  with  experimental  results.  It  appears  that 
the  detrimental  influence  of  the  moisture  is  the 
degradation  of  the  matrix  shear  modulus  at  high 
compressive  stresses  with  high  moisture  concentrations. 
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BACKGROUND 


INTRODUCTION 

Fiber  reinforced  composite  materials  have  become 
increasingly  popular,  in  part,  due  to  their  high  strength 
to  weight  ratios  and  the  ability  to  orient  the 
reinforcement  in  the  principal  stress  directions.  The 
primary  applications  of  such  materials  have  been  in  the 
aircraft  and  sporting  industries  with  typical  products 
usually  involving  thin  sections  with  tensile  loadings. 

New  applications  for  these  materials  are  currently  being 
explored  for  naval  and  biomedical  hardware  where  the 
products  involve  thicker  sections  loaded  in  compression 
while  exposed  to  a  moist  environment. 

The  usefulness  of  the  composite  materials  database  as 
it  currently  exists  for  these  new  applications  is  unclear 
because  studies  conducted  over  the  past  twenty  years  have 
focused  primarily  on  the  behavior  of  thin  sections  loaded 
in  tension.  Composite  materials  often  do  not  exhibit  the 
same  behavior  in  compression  as  in  tension  [1,2]. 
Differences  in  elastic  moduli,  strengths  and  failure  modes 
have  been  observed.  Fabrication  quality  also  becomes  more 
important  if  the  component  is  to  be  loaded  in  compression. 
Compressive  strengths  are  very  sensitive  to  composite  void 
contents  and  initial  fiber  waviness  [3,4, ',6,7],  and 
void  contents  of  less  than  1%  to  2%  are  required  to 
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avoid  significant  reductions  in  compression  strength  [8], 

Current  compression  test  methods  for  composites  have 
been  developed  for  thin  sections  based  on  the 
requirements  of  the  aircraft  industry.  These  methods 
transfer  the  load  by  shear  from  the  grips  to  the 
specimen.  Characterization  of  the  compression  response 
of  thick  section  composites  using  these  methods  is 
unlikely.  The  very  large  loads  required  for  failure  of 
the  thick  sections  eliminate  the  option  of  transferring 
the  load  by  shear  from  the  grips  to  the  specimen.  New 
tet.*-  methods  may  have  to  be  developed  to  assess  the 
effect  of  thickness  on  compression  behavior  [9], 

Environmental  conditioning  introduces  additional 
complications.  The  typical  conditioning  procedure  used  by 
the  aircraft  industry  requires  full  saturation  of  the 
specimen  in  a  hot,  humid  environment  [10,11,12].  However, 
this  procedure  may  be  unrealistic  for  thick  sections  with 
applications  in  environments  that  require  total  immersion 
in  seawater  at  or  near  room  temperature.  Full  saturation 
of  thick  sections  at  these  low  temperatures  is  not  feasible 
due  to  the  extremely  long  period  of  time  required  for 
saturation . 

Fabrication,  environmental  conditioning  and 
compression  testing  of  specimens  with  thin  sections  is 
easier  and  less  expensive  than  comparable  procedures  for 
thick  laminates.  Nevertheless,  if  significant  differences 
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are  observed  between  the  mode  of  failure  for  thin  and 
thick  specimens,  then  it  is  clear  that  thin  specimens 
cannot  be  used  to  characterize  properties  of  thick  section 
laminates . 

The  first  objective  of  this  study  is  to  address  these 
issues  by  examining  the  changes  in  the  failure  mode  and/or 
strength  between  thin  1.6  mm  (0.0625  in.)  and  thick 
6.35  mm  (0.25  in.)  graphite/epoxy  laminates  loaded  in 
uniaxial  compression.  The  second  objective  is  to  compare 
the  fracture  behavior  for  both  dry  and  "wet"  (moisture 
conditioned)  specimens. 
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LITERATURE  REVIEW 

COMPRESSION  TESTING  OF  COMPOSITE  MATERIALS 

Compression  testing  of  composite  materials  usually 
follows  the  procedures  outlined  in  ASTM  D3410-87  [13]. 

Two  fixtures  recommended  in  this  standard  are  the  Celanese 
and  IITRI  (Illinois  Institute  of  Technology  Research 
Institute)  compression  test  fixtures  [14,15].  These  two 
fixtures  are  illustrated  in  Figures  1  and  2.  With  both 
fixtures,  the  compressive  load  is  transmitted  by  shear  from 
the  wedged  grips  to  the  specimen.  The  IITRI  fixture 
differs  from  the  Celanese  primarily  in  the  design  of  the 
wedged  grips  [13].  Trapezoidal  wedges  are  used  in  the 
IITRI  fixture  which  maintain  surface  contact  between  the 
grips  and  the  alignment  blocks  for  specimens  with  tabbed 
end  thicknesses  of  1.5  mm  (0.06  in.)  to  12.7  mm  (0.50  in.) 
[14].  The  conically  shaped  wedges  of  the  Celanese  fixture 
limit  the  tabbed  end  thickness  to  a  maximum  of  4.06  mm 
(0.16  in.)  because  greater  thicknesses  cause  a  line  instead 
of  a  surface  contact  between  the  conical  wedge  and  tapered 
collet  [ 14 ] . 

Recommended  specimen  geometries  for  either  fixture  are 
also  given  in  ASTM  D3410-87,  and  these  dimensions  depend 
upon  the  material  being  tested.  For  example,  the 
thicknesses  of  a  unidirectional  laminate  of  graphite/epoxy 
should  range  between  1.5  to  3.0  mm  (0.06  to  0.120  in.). 
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Specimen  widths  and  gage  lengths  are  specified  at  6.35  mm 
(0.25  in.)  and  12.70  mm  (0.50  in.),  respectively.  Tabs  are 
required  for  all  specimens  and  should  be  63.5  mm  (2.50  in.) 
long  and  6.35  mm  (0.25  in.)  wide  with  a  6.35  mm  (0.25  in.) 
taper  toward  the  specimen  gage  area.  Design  and  precision 
in  machining  the  specimen  used  with  this  test  method  is 
critical.  To  minimize  bending  strains  due  to  specimen 
eccentricity  requires  that  the  tab  surfaces  be  ground  flat 
and  parallel  to  within  0.05  mm  (0.002  in.)  [16].  In 
addition,  the  specimen  gage  lengths  recommended  in  / 

D3410  must  be  reviewed  for  new  composite  materials  with 
varying  degrees  of  anisotropy.  The  gage  length  must  be 
long  enough  to  accomodate  the  stress  concentration  effects 
introduced  at  the  tab  ends,  but  short  enough  to  avoid  Euler 
buckling.  The  specimen  length  necessary  to  attenuate 
stress  concentrations  due  to  the  specimen  ends  and/or  grips 
for  anisotropic  materials  can  exceed  several  specimen 
thicknesses  [17]. 

Composite  materials  have  also  been  tested  in 
compression  using  the  sandwich  beam  specimen  shown  in 
Figure  3  [18].  This  specimen  consists  of  two  skins  of 
composite  material  bonded  to  a  honeycomb  core.  The  skin 
that  is  exposed  to  compression  when  bending  loads  are 
applied  to  the  specimen  is  made  thinner  than  the  tension 
skin  to  force  failures  in  the  compression  layer. 

Each  of  the  preceding  methods  have  disadvantages.  The 
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requirement  for  specimen  tabs  significantly  complicates 
specimen  preparation  for  the  IITRI/Celanese  methods.  In 
addition,  testing  of  high  strength  materials  or  thicker 
cross-section  specimens  is  limited  by  the  shear  strengths 
of  available  adhesives  which  are  used  to  bond  the  tabs  to 
the  specimens  [14],  Typical  adhesive  shear  strengths 
range  from  35  to  42  MPa  (5  to  6  ksi)  limiting  compressive 
loads  to  33.4  kN  (7500  lbs.)  for  a  63.5  mm  (2.50  in.)  by 
6.35  mm  (0.25  in.)  attachment  area  [19].  The  specimen 
length  required  for  the  attenuation  of  stress 
concentrations  due  to  the  specimen  ends  and/or  grips 
also  increases  with  thickness  [17].  This  requirement 
complicates  specimen  design  as  analyses  have  to  be 
performed  for  each  specimen  thickness  to  determine  the 
appropriate  decay  length.  Finally,  the  sandwich  beam 
specimen  is  applicable  only  to  very  thin  laminates  1.020 
mm  (0.040  in.),  and  requires  a  large  amount  of  composite 
material  for  fabrication  [14,  18]. 

Loading  of  composites  in  compression  directly  through 
the  specimen  ends  is  advantageous  as  the  requirement  for 
end  tabs  and  limits  on  cross-sectional  areas  are 
eliminated.  Another  standard,  ASTM  D695,  a  compression 
test  for  rigid  plastics,  recommends  loading  the  ends  of  a 
"dogbone"  shaped  specimen  which  is  bolted  in  a  fixture 
that  supports  the  specimen  length.  This  fixture  and 
specimen  are  shown  in  Figures  4  and  5  [20].  The  specimen 
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configuration  is  only  for  materials  with  thicknesses  of 
3.2  mm  (0.125  in.)  or  less,  and  has  not  been  found 
acceptable  for  the  high  modulus  fiber  reinforced  composite 
materials  considered  in  this  investigation  [15].  Examples 
of  compression  fixtures  for  thicker  specimens  are  also 
provided  as  guidelines  in  this  test  standard.  To  date, 
fixtures  which  apply  loads  through  the  specimen  ends 
appear  to  be  independently  developed  [9,14,21,22].  A 
distinguishing  feature  of  these  designs  is  that  some  type 
of  lateral  support  is  usually  required  to  prevent  end 
brooming,  particularly  in  unidirectional  laminates  [14]. 
For  example,  Barker  and  Balasundaram  [21]  developed  a 
fixture  consisting  of  two  steel  clamping  blocks  which 
supported  all  four  sides  of  the  specimen  ends  over  a 
length  of  40  mm  (1.57  in.).  The  specimen  ends  were 
ground  flat  and  parallel,  and  the  compression  load  was 
applied  directly  to  these  ends.  In  using  this  fixture, 
they  observed  that  approximately  20%  of  the  total  load 
had  to  be  transmitted  by  shear  forces  developed  by  the 
steel  blocks.  If  all  the  load  was  transmitted  through 
the  ends,  premature  failures  occurred  within  the  clamping 
blocks.  If  the  shear  forces  exceeded  20%  of  the  load, 
premature  failures  occurred  at  the  edge  of  the  clamp. 

A  similar  fixture  was  designed  by  Camponeschi  [9]  at 
the  David  Taylor  Research  Center.  Two  steel  clamping 
blocks  support  the  specimen  at  each  end  as  illustrated  in 
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Figure  6.  The  bolt  torque  required  to  prevent  end 
brooming  due  to  Poisson  effects  is  dependent  on  the 
specimen  thickness,  and  must  be  calculated.  Also,  tabs 
were  adhesively  bonded  to  each  specimen  over  the  clamping 
area. 

Comparisons  of  test  results  obtained  using  the  three 
preceding  test  methods,  IITRI/Celanese,  sandwich,  and  end 
loading,  show  that  the  longitudinal  elastic  modulus  is 
independent  of  the  compression  test  method  used  [14,15]. 
The  results  for  compressive  strength  determined  with  the 
IITRI/Celanese  and  sandwich  beam  methods  are  comparable  if 
the  specimens  are  thin  and  properly  prepared  [15].  Some 
difficulties  have  been  reported  using  the  Celanese  method 
where  slippage  occurred  between  the  grips  and  the  tabs  on 
graphite/epoxy  specimens.  However,  slippage  was  not 
observed  with  the  IITRI  fixture.  The  problem  was 
attributed  to  the  shallow  serrations  of  the  Celanese 
fixture  as  compared  to  the  coarser  serrations  found  in  the 
IITRI  grips  [14]. 

Variations  in  the  compression  strength  and  the  failure 
mode  obtained  with  the  IITRI  fixture  have  been  noted  with 
changes  in  the  gage  length  to  thickness  ratio  [23],  Smoot 
reported  that  the  highest  compressive  strengths  for 
unidirectional  graphite/epoxy  occurred  with  a  length  to 
thickness  ratio  of  8.7.  The  failure  mode  corresponding  to 
that  strength  was  at  a  transition  between  angular  shear 
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across  the  specimen  width  and  brooming  within  the  gage 
section.  Schematics  of  these  failure  modes  are  shown  in 
Figure  7  [23].  Exceeding  a  length  to  thickness  ratio  of 
8.7  resulted  in  a  rapid  decrease  in  failure  strength, 
triggered  by  specimen  Euler  buckling.  Strengths  also 
dropped,  but  less  rapidly,  if  a  length  to  thickness  ratio 
less  than  8  was  used.  It  should  be  noted  that  the  ASTM 
Standard  D3410  recomends  a  length  to  thickness  ratio  in  the 
range  of  4.2  to  8.3. 

Test  methods  using  end  loading  have  met  with  limited 
success.  Lower  strengths  for  unidirectional  laminates 
have  been  reported  using  these  methods  as  compared  to  the 
other  compressive  test  methods  [9,15,22].  The  lower 
strengths  were  attributed  to  premature  failures  by  end 
brooming.  This  type  of  failure  may  be  due  to  Poisson 
effects,  or  to  stress  concentrations  caused  by  uneven 
fiber  ends  or  transverse  friction  loads  on  the  ends  of  the 
specimen  [14,15].  However,  Barker  and  Balasundaram,  using 
the  end  loading  fixture  previously  described,  obtained 
compressive  strengths  of  1470  MPa  (210  ksi)  for 
unidirectional  graphite/ epoxy .  This  correlates  very  well 
with  strengths  obtained  by  the  other  methods  [15,23]. 

The  use  of  end  loading  techniques  with  cross  ply  or 
angle  ply  laminates  has  been  more  successful  [9,22]. 
Strengths  of  unidirectional  laminates  of  graphite/epoxy 
estimated  from  cross  ply  compression  test  results  (6.35  mm 
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thick  specimens)  using  Camponeschi 's  fixture  agreed  well 
with  those  determined  using  the  IITRI/Celanese  methods 
[9].  These  more  satisfactory  results  may  be  due  to  the 
transverse  support  provided  by  the  90°  reinforcement  in 
the  cross  ply  laminate  [24].  However,  as  the  specimen 
thicknesses  were  increased  from  6.35  mm  (0.25  in.)  to  25.4 
mm  (1.00  in.),  the  cross  ply  laminate  strengths  dropped  by 
21%.  Since  laminates  of  these  thicknesses  have  rarely 
been  tested,  it  is  not  possible  to  conclude  if  the  loss  of 
strength  is  due  to  material  changes,  or  to  the  test 
method . 

Pressure  cylinders  fabricated  with  composite  materials 
with  wall  thicknesses  exceeding  12.7  mm  (0.50  in.)  provide 
one  source  of  comparison.  The  wall  stresses  in  thick 
walled  graphite/epoxy  cylinders  at  failure  were 
significantly  below  typical  flat  laminate  strengths 
[9  25].  This  trend  was  not  observed  in  fiberglass 
rexnforced  cylinders  [26,27]. 

For  a  more  complete  review  of  compression  test  methods 
for  composite  materials,  the  reader  is  referred  to 
Reference  28. 
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COMPRESSION  FAILURE  THEORIES 

Since  the  1960's,  numerous  researchers  have  attempted 
to  analytically  describe  and  quantify  the  failure  of 
composite  materials  under  compressive  loads.  Development 
of  compression  failure  theories  has  been  difficult  due  to 
the  multiplicity  of  failure  modes  seen  with  variations  in 
constituent  properties  and  laminate  construction  [1].  If 
these  failure  theories  are  grouped  according  to  the 
particular  failure  mode  addressed,  three  general 
categories  of  failure  theories  emerge  —  fiber 
microbuckling,  kink  band  formation,  and  combined 
failure  modes.  A  brief  characterization  of  typical 
compression  failure  theories  from  each  of  the  three 
categories  is  given  below. 

Historically,  the  failure  mode  most  frequently 
applied  to  unidirectional  composites  subjected  to 
compressive  stress  is  that  of  fiber  microbuckling.  Rosen 
[29]  proposed  one  of  the  earliest  theories  based  on  fiber 
microbuckling  where  the  fibers  and  matrix  were  modeled  as 
Timoshenko  columns  supported  by  an  elastic  foundation. 
Expressions  for  the  compression  strength  were  derived  for 
both  the  extensional  and  shear  buckling  modes.  In  the 
extensional  mode,  the  fibers  buckle  out  of  phase,  and  in 
the  shear  mode,  they  buckle  in  phase.  Equation  (1)  shows 
Rosen's  expression  for  the  shear  mode  of  failure  in 
compression,  which  was  the  mode  most  commonly  observed  in 


unidirectional  coroosites. 

°c=glt  (D 

oc=  compression  strength 

Gu=  inplane  shear  modulus  of  composite 

Compression  strength  values  obtained  with  this 
theoretical  expression  significantly  overpredicted 
experimental  results. 

Additional  fiber  microbuckling  theories  including 
three  dimensional  models,  initial  ply  waviness  and 
material  nonlinearities  were  developed  in  an  attempt  to 
improve  correlations  between  theory  and  experimental 
observations  [30,31,32].  Using  a  3-dimensional  model  and 
energy  methods,  Greszczuk  derived  two  expressions  for  the 
microbuckling  stress  of  a  composite  based  on  the  fixity  of 
the  fiber  ends. 


oc=Glt+( — -f-1)  (y)2  Simply  supported 

(2) 

ac=GLT+7T2E(Vf  (y) 2  Fixed 

(3) 

Glt=  in-plane  shear  modulus  of  composite 
E,=  longitudinal  modulus  of  fiber 
Vf=  fiber  volume  fraction 
r=  fiber  radius 
1=  specimen  length 

It  was  noted  that  these  expressions  do  not  account  for 
nonlinear  matrix  behavior  or  the  initial  fiber 
deformations  which  could  significantly  affect  the 
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microbuckling  loads.  Good  agreement  was  demonstrated 
between  this  theory  and  experiments  with  composites 
fabricated  from  low  modulus  resins.  However,  for 
composites  with  high  modulus  resins,  the  theory  is  not 
applicable  because  failure  is  not  due  to  microbuckling . 
This  fact  led  to  the  development  of  another  relation  for 
the  compression  strength  based  on  a  rule  of  mixtures 
formulation  [30]. 

°c=<W-M] L"v()  (4> 

of=  fiber  compressive  strength 

or=  matrix  compressive  stress  corresponding 

o, 

to  a  strain  of  W- 

The  effects  of  initial  fiber  waviness  and  matrix 
nonlinearities  were  addressed  by  Davis  who  modeled  fiber 
and  matrix  layers  as  a  multilayered  Timoshenko  beam  [33]. 
Critical  strength  expressions  for  failure  by  delamination 
and  failure  by  shear  instability  were  developed  and  these 
relations  correlated  well  with  experimental  results  for 
boron/epoxy  laminates.  The  critical  failure  load  for 
shear  instability  was  a  function  of  the  shear  moduli  and 
volume  fractions  of  both  the  fiber  and  matrix  as  shown 
below. 

GfGm 

°C  VfGm+VmG(  ^ 

Gf=  shear  modulus  of  fiber 
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Gm=  shear  modulus  of  matrix 
Vm=  matrix  volume  fraction 


However,  Davis  cautions  that  the  matrix  shear  modulus 
used  in  Equation  (5)  must  correspond  to  the  shear  stress- 
strain  state  in  the  matrix,  so  that  the  nonlinearity  of 
the  shear  modulus  of  the  matrix  is  reflected  in  the 
solution.  The  stress-strain  state  in  the  matrix  is 
determined  from  a  finite  difference  solution  for  the 
multilayered  Timoshenko  beam  subjected  to  an  axial  load. 

Hahn  [1]  also  developed  an  expression  for  the  critical 
microbuckling  stress  of  a  composite  incorporating  both 
initial  fiber  waviness  and  matrix  nonlinearities. 


°c~vf  „Efdf  7i f 

1+htri-r> 


(6) 


ob=  fiber  flexural  strength 
df=  fiber  diameter 

~Y=  amplitude  to  length  ratio  of 

initial  fiber  imperfection 


Experimental  validation  of  Equation  (6)  showed  reasonable 
agreement  for  composites  with  resin  tensile  moduli  between 
3  and  5  GPa  (.429E6  psi  and  .716E6  psi)  . 

The  second  group  of  failure  theories  which  are  closely 
related  to  fiber  microbuckling,  focus  on  the  development 
of  kink  bands  in  unidirectional  composites  under 
compressive  loads  (See  Figure  8).  Hahn  [1]  found  that 
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fiber  microbuckling  occurred  in  composites  with  low 
modulus  resins,  while  fiber  kinking  was  the  typical 
failure  mode  for  composites  with  high  modulus  resins. 
Microscopically,  fiber  microbuckling  was  characterized  by 
high  bending  strains  and  the  consequent  development  of 
multiple  fiber  breaks.  Hahn  defined  failure  in  the 
kinking  mode  with  two  breaks  per  fiber,  where  the  short 
broken  segments  align  in  a  parallel  pattern  rotated  at  an 
angle  less  than  90°  to  the  applied  compressive  load. 
Typical  kink  band  angles  for  graphite/ epoxy  ranged  from 
22°  to  26°.  No  attempt  was  made  to  develop  a  theoretical 
expression  for  the  compression  strengths  based  on  the  kink 
band  failure  mode. 

Kink  band  formation  under  compressive  load  was  also 
described  by  Budiansky  [34],  and  incorporated  the  effects 
of  matrix  plasticity  and  initial  fiber  waviness.  The 
compression  strength  based  on  this  theory  was  a  function 
of  the  kink  band  angle,  the  effective  transverse  shear 
modulus  of  the  composite  and  the  transverse  Young's 
modulus  as  shown  in  Equation  (7) . 

o=  C^]  (GLT+ETtan(3)  (7) 

y  =  strain  corresponding  to  shear  yield  stress 

7i=  initial  imperfection 

Glt=  inplane  shear  modulus  of  composite 

Et=  transverse  modulus  of  composite 

[5=  kink  band  angle 
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From  the  preceding  discussion,  it  is  clear  that  a 
single  expression  for  the  compression  strength  of  a 
composite  does  not  predict  the  multiple  failure  modes 
observed  under  compressive  loads.  This  fact  led  a  number 
of  researchers  to  develop  multiple  mode  failure  theories. 
Sinclair  and  Chamis  [35]  developed  four  expressions  for 
compressive  strengths  each  linked  to  a  different  failure 
mode.  Their  models  addressed  fiber  microbuckling, 
delamination,  bending,  and  compression  failure  of  the 
fibers.  A  similar  approach  was  used  by  Piggott  [36]  who 
derived  six  different  failure  models. 

More  recently,  Shuart  [37]  developed  a  method 
considering  three  mechanisms  leading  to  failure  under 
compressive  loads.  These  included  matrix  compression, 
interlaminar  shearing  due  to  out  of  plane  waviness  in  the 
ply,  and  inplane  matrix  shearing  which  was  affected  by 
inplane  fiber  waviness  and  fiber  scissoring.  The 
uniqueness  of  Shuart 's  approach  is  that  it  addressed 
compression  failures  in  both  unidirectional  and  angle  ply 
laminates.  The  composite  laminate  was  modeled  as  a  series 
of  Kirchhoff  plates  and  elastic  foundations  with  initial 
out  of  plane  imperfections  having  the  same  shape  as  the 
critical  short  wavelength  buckling  mode.  Nonlinear 
equilibrium  equations  for  the  laminate  were  derived  using 
von  Karman  nonlinear  strain-displacement  relations  and  the 
principle  of  stationary  potential  energy.  Solution  of 
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these  equations  assuming  linear  elastic  material 
properties  yielded  displacements,  strains,  and  stresses. 

Experimental  investigations  conducted  by  Shuart  showed 
that  one  failure  mechanism  predominated  for  a  particular 
range  of  angles  in  the  angle  ply  laminates.  For  example, 
interlaminar  shearing  was  the  dominant  failure  mode  when 
the  longitudinal  axes  of  the  plies  were  oriented  from  0° 
to  15°  from  the  loading  axis.  Theoretical  predictions  of 
failure  mode  and  strength  generally  correlated  well 
(within  5%)  with  experimental  observations  for  laminates 
with  plies  oriented  between  45°  and  90°.  However,  for 
unidirectional  laminates,  the  failure  strengths  predicted 
were  significantly  less  than  the  experimentally  determined 
strengths. 

The  compression  failure  theory  of  Davis  [6]  was 
chosen  for  this  study  to  compare  the  difference  in 
compression  strengths  between  wet  and  dry  specimens  for 
several  reasons.  First,  is  that  the  author  demonstrated 
good  agreement  between  theory  and  experiment  for 
unidirectional  boron/epoxy  laminates  which  failed  by  shear 
instability  (as  defined  in  Reference  [6]).  This  is 
important  because  the  response  of  the  (02/90)s  laminates 
in  this  study  will  be  dominated  by  the  0°  plies,  and  shear 
instability  has  been  an  observed  failure  mode  for  these 
laminates  in  previous  studies  [9,23].  Shuart  [22]  also 
predicted  interlaminar  shearing  as  a  failure  mode  for 
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unidirectional  laminates,  but  his  predicted  compression 
strengths  were  significantly  lower  than  those  commonly 
observed  for  unidirectional  AS4/3 501-6  laminates. 

Secondly,  the  form  of  the  failure  theory,  as  shown  in 
Equation  (5) ,  is  representative  of  many  of  the 
microbuckling  based  theories  and  Budiansky's  [34]  kink 
band  theory.  For  example,  if  the  shear  modulus  (G^)  of  a 
multilayer  beam  is  derived  by  considering  the  shear 
modulus  of  each  layer  as  springs  in  parallel,  then 
Equation  (5)  states  that  the  compression  strength  of  a 
composite  is  equivalent  to  its  shear  modulus.  This 
conclusion  is  the  same  as  Rosen's  [29],  and  a  lower  bound 
for  Greszczuk's  [30]  microbuckling  theories. 

And  lastly,  since  the  shear  modulus  of  a  composite  is 
highly  nonlinear,  some  methodology  must  be  proposed  to 
choose  a  modulus  for  use  in  the  failure  theory.  Davis  [6] 
recommended  such  an  approach  including  the  effects  of 
nonlinear  matrix  behavior  and  fiber  waviness.  Davis  based 
his  approach  on  Timoshenko  beam  theory  and  solved  the 
resultant  equations  with  finite  difference  techniques,  but 
other  numerical  techniques  such  as  finite  element  methods 
would  be  equally  applicable  to  this  problem. 
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DIFFUSION  BEHAVIOR  IN  COMPOSITE  MATERIALS 


Understanding  the  effects  of  moisture  on  composite 
mechanical  behavior  requires  the  ability  to  estimate  the 
amount  of  water  a  composite  will  absorb  under  a  specified 
set  of  environmental  conditions.  Most  of  the  data 
currently  available  regarding  composite  materials  and 
moisture  absorption  has  been  generated  by  the  aircraft 
industry.  Graphite/epoxy  is  the  material  most  frequently 
cited,  and  the  environmental  exposure  is  typically  hot  and 
humid  [10,11,38].  Specimens  are  usually  exposed  to  this 
environment  until  full  saturation  occurs  which  is 
manifested  by  no  further  increases  in  specimen  weight  with 
continuing  exposure  [38].  Much  less  information  is 
available  on  the  effects  of  immersion  as  compared  to 
humidity  conditioning,  and  moisture  gradients  through  the 
thickness  as  compared  to  a  constant  moisture 
concentration.  These  issues  are  pertinent  to  new 
applications  for  composite  materials.  Moisture  diffusion 
in  composites  is  further  complicated  by  the  number  of 
parameters  that  can  affect  the  process,  which  include  void 
content,  fiber  orientation,  stress  state,  resin  material, 
and  degree  of  cure  [39,40].  If  it  is  assumed  that  the 
fiber  is  impervious,  all  the  moisture  absorbed  by  a 
composite  is  attributed  to  the  resin  material  [41].  The 
rate  of  diffusion  of  water  into  a  polymer  material  is  then 
modeled  by  Fick's  Law  [42]: 
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3C_n82C 

0Z2 


(8) 


C=moisture  concentration 

D=diffusion  coefficient  through  thickness 
z=thickness  coordinate 


The  solution  of  this  partial  differential  equation  for  a 
system  with  an  initial  uniform  moisture  concentration  of 
Cq/  and  moisture  concentrations  of  C1  and  C at  the  -h/2 
and  +h/2  (see  Figure  9)  surfaces  respectively  is  [43,44] 


7  o  ®  C,cos(nn)  -C.  . 

C  ( z  ,  t)  =C,+  (C2-C,) §+f  2— - n - sin  (^)  e h 

“  rul  11 

,4C. ?  1  - , n(2m+l)»z 

*  i„2m+1  h 


2,2 


(9) 


i  or  experimental  purposes, 
the  per  cent  weight  gain. 


it  is  helpful  to  calculate  M, 


(10) 


Ww=  weight  of  wet  material 
WD=  weight  of  dry  material 


which  can  also  be  expressed  as  [38] 

h 

M=G  Ji^caz 

2 


(11) 


G=  time  dependent  parameter 
defined  below 


The  integration  of  Equation  11  gives  [45] 
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M-M. 


-(2j*nV(2|) 

!.Iye _ 

*2h  (2 j+i) 2 


(12) 


Mt=initial  moisture  content 
M^maximum  moisture  content 

D=diffusion  coefficient  through  the  thickness 

t=time 

h=thickness 


Springer  [38]  approximates  this  solution  for  a  material 
exposed  to  moisture  on  two  sides  by 

G=l-exp[-7.3(^)'75]  (13) 

h 

A  composite  material  formed  of  water  absorbing  resin 
and  an  impervious  fiber  will  absorb  less  moisture  than  the 
pure  resin-  The  diffusion  coefficient  of  moisture  through 
the  thickness  will  be  a  function  of  fiber  angle  and  fiber 
volume  fraction.  An  expression  for  this  coefficient  has 
been  developed  by  Shen  and  Springer  [38] 

D=D„cos2a+D22sin2a  (14) 


D,rDr(l-Vf) 

D22=Dr(l-2^) 


Dn=dif fusivity  parallel  to  fibers 
D22=dif fusivity  perpendicular  to  fibers 
Dr=resin  diffusion  coefficient 
a=orientation  of  fibers  with  respect 
to  z  axis, 90°  (See  Figure  10) 

Vf=  fiber  volume  fraction 
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If  a  composite  specimen  is  exposed  to  moisture  without 
coating  its  edges,  correction  for  the  edge  effects  must  be 
incorporated  into  the  calculation  of  the  diffusion 
coefficient.  For  a  unidirectional  composite  oriented  as 
shown  in  Figure  10: 


D2=Dncos2a  1  L'22sin2a 

(15) 

Dy=Ducos2p+D22sin2p 

(16) 

Dx=Dncos2(  90-P)  +D22sin2  ( 90-P) 

(17) 

p=orientation  of  fibers  with  respect 
to  y  axis 


In  the  early  stages  of  diffusion,  i.e.  when  the  ratio  of 
the  moisture  absorbed  to  the  maximum  moisture  content 
varies  linearly  with  the  square  root  of  the  time  of 
exposure,  there  is  no  significant  interaction  between  the 
diffusion  of  the  different  sides.  Equation  (8)  can  then  be 
applied  independently  to  each  edge  of  the  composite 
resulting  in  a  diffusion  coefficient  corrected  for  edge 
effects  [ 38 ] . 

h  [Dl  h  ID*  ? 

d=dz(1+t^t^)  <18> 


For  a  laminated  composite,  Equations  15  to  17  become 


D2=D22-D,(1-2\|5) 

Ih  cos2P  Ih,sin2P 
Dy=D»  ZhJ  ^"D22  ihj  ’ 

Ih,sin2P  Ih,cos2P 

p.-D>.  ih,  1+p^ — j 


(19) 

(20) 

(21) 
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hj=thickness  of  jlh  ply 

(l=angle  of  jth  ply  with  respect  to  y  axis 

The  preceding  equations  are  based  on  Fickian 
diffusion.  Non-Fickian  moisture  absorption  behavior  may 
be  observed  in  several  instances  such  as  when  moisture 
exposure  occurs  at  elevated  temperatures  and  microcracking 
develops  [46,47]  or  when  exposure  occurs  while  the 
material  is  being  mechanically  stressed  [40].  Some 
researchers  have  reported  non-Fickian  diffusion  behavior 
for  the  3501-6  resin  following  exposure  above  certain 
temperatures  and  humidities  [46,47].  Collings  and  Copley 
found  non-Fickian  diffusion  for  conditioning  temperatures 
above  60°  C  (140°  F)  or  humidities  above  96%  RH  [52]. 
However,  Shirrell  and  Sandow  observed  non-Fickian 
diffusion  in  the  3501  resin  only  at  combined  conditions  of 
elevated  temperature  (82°  C)  and  humidity.  All  the 
specimens  that  exhibited  this  anomaly  had  surface  and  edge 
microcracks  [46].  In  such  cases,  the  water  absorption  may 
be  higher  than  calculated  by  the  Fickian  model.  An 
increase  in  the  fraction  of  voids  will  also  increase  the 
moisture  weight  gain. 

Matrix  plasticization  is  only  one  effect  of  composite 
moisture  absorption.  Matrix  swelling  with  moisture 
can  also  be  significant.  The  swelling  strains  due  to 
moisture  are  treated  analogously  to  thermal  strains. 
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em=Pm(CrC() 


(22) 


em=rnoisture  swelling  strain 
Pm=moisture  expansion  coefficient 
Cf=final  moisture  concentration 
C,=initial  moisture  concentration 

Differential  stresses  can  develop  in  a  composite  with  a 
moisture  gradient  through  the  thickness  as  the  outer 
layers  expand  and  pull  on  the  inside  layers  [48,49].  The 
effect  of  these  swelling  stresses  may  be  positive  or 
negative  depending  on  whether  they  relieve  or  add  to  the 
residual  stresses  in  the  laminate. 
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CONDITIONED  RESULTS 


The  available  data  on  mechanical  property  changes  in 
composite  materials  due  to  environmental  conditioning  has 
been  generated  primarily  by  the  aircraft  industry. 
Consequently,  the  environmental  conditioning  reflects 
aircraft  operating  conditions  i.e.  thin  composite  sections 
exposed  to  hot,  humid  environments. 

Differences  in  exposure  method,  moisture  content, 
materials,  laminate  orientations,  temperature  and  loading 
condition  make  it  difficult  to  draw  general  conclusions 
regarding  the  effects  of  environmental  conditioning  on 
composite  mechanical  properties.  Decreases  in  compression 
strength  ranging  from  11%  to  20%  have  been  reported  for 
thin  graphite/ epoxy  laminates  humidity  conditioned  at  room 
temperature  to  50%  maximum  saturation  and  higher  [21,50,51 
52].  The  longitudinal  elastic  modulus  of  unidirectional 
laminates  appears  to  remain  unchanged  with  moisture  at 
room  temperature  [38,50],  but  the  transverse  modulus  drops 
by  12%  [38]  . 

Failure  mode  differences  between  dry  and  wet  specimens 
are  also  of  interest  as  they  may  provide  insight  into  the 
differences  in  compression  strengths  between  these 
specimens.  Grimes  and  Adams  [50]  found  a  20%  drop  in  the 
static  compression  strength,  and  observed  similar  gross 
failure  modes  in  the  dry  and  wet  (unidirectional 
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graphite/epoxy  laminates  with  a  1%  moisture  content  by 
weight) .  However,  the  wet  specimens  showed  slightly  more 
delamination  than  the  dry  specimens.  Their  specimens  were 
tested  with  the  IITRI  compression  fixture.  Barker  and 
Balasundaram  [21]  tested  2.00  mm  (0.079  in.)  thick 
laminates  by  end  loading  in  compression,  and  found  changes 
in  failure  mode  in  the  0°  lamina  of  a  [(90/0)4]s  laminate 
from  shear  in  the  dry  material  to  fiber  microbuckling  in 
the  wet  material.  These  changes  in  failure  mode  were 
accompanied  by  a  12.5%  drop  in  compression  strength. 
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EXPERIMENTAL  PROCEDURE 


SPECIMEN  FABRICATION  AND  PREPARATION 

Four  graphite/epoxy  panels  were  fabricated  from 
AS4/3501-6  prepreg  (.005"  nominal  thickness).  The  fiber 
and  epoxy  properties  are  shown  in  Tables  1  and  2.  Two 
panels  were  prepared  with  48  plies  to  give  a  nominal  6.35 
mm  (0.25  in.)  thickness  when  cured.  The  remaining  two 
panels  had  12  plies  and  a  1.6  mm  (0.0625  in.)  cured 
thickness.  The  time-temperature  profile  used  in  curing  is 
shown  in  Figure  11.  A  symmetric  0^/90  layup  was  selected 
instead  of  a  unidirectional  layup  to  avoid  difficulties 
due  to  transverse  brooming  expected  with  testing  thick 
unidirectional  coupons.  The  90°  reinforcement  provides 
transverse  constraint  which  aids  development  of  the  full 
material  strength. 

Each  panel  was  inspected  for  delaminations  and  voids 
with  an  ultrasonic  C-scan.  Tag  end  samples  were  used  for 
acid  digestion  to  determine  fiber  volume  fractions  of  62% 
and  66%  for  the  thick  and  thin  panels,  respectively  [53]. 
Void  contents  using  this  technique  were  typically  less  than 
1%.  Specimens  with  the  dimensions  shown  in  Figure  12 
were  cut  from  each  panel.  The  ends  of  each  specimen  were 
ground  flat  and  parallel  to  a  tolerance  of  +0.0254  mm 
(+0.001  in.)  to  minimize  difficulties  in  loading  due  to 
uneven  fiber  ends.  The  gage  length-to-thickness  ratio  (l:t 
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ratio)  for  the  thick  specimens  was  3:1.  This  1/t  provides 
a  safety  factor  of  2.0  against  Euler  buckling  assuming 
clamped  ends  and  including  transverse  shear  effects  (see 
Appendix  A) .  The  same  1/t  ratio  for  the  thin  specimens 
means  the  gage  length  would  be  4.7625  mm  (0.1875  in.). 

This  1/t  ratio  was  used  initially,  but  resulted  in  low 
strength  failures  at  the  specimen's  exit  from  the  grips. 
Consequently,  the  1/t  ratio  for  the  thin  specimens  was 
increased  to  6.4:1  so  that  stress  concentration  effects 
from  the  grips  would  not  extend  throughout  the  entire  gage 
length.  The  buckling  safety  factor  for  this  1/t  ratio  was 
1.1.  Selected  specimens  of  each  group  were  strain  gaged 
for  modulus  measurements  with  a  single  gage  or  back  to 
back  gages  to  monitor  specimen  alignment.  CEA-06-125UW- 
350  gages  were  used  on  the  6.35  mm  (0.25  in.)  thick 
specimens  with  a  gage  length  of  19.05  mm  (0.75  in.).  The 
strain  gages  were  applied  immediately  prior  to  testing 
with  a  cyanoacrolate  adhesive.  The  same  procedure  was 
followed  for  the  thin  specimens  using  CEA-06-062UT-350 
gages  to  accomodate  the  shorter  gage  length. 

One  objective  of  this  investigation  was  to  compare  dry 
and  wet  specimen  behavior.  This  required  moisture 
conditioning  of  half  the  specimens.  The  specimens  to  be 
moisture  conditioned  were  first  dried  repeatedly  in  a  50°  C 
(122°  F)  oven  until  there  was  less  than  a  0.5%  change  in 
specimen  weight  before  and  after  drying.  After 
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conditioning,  the  specimens  were  stored  in  room  temperature 
^ater  until  testing.  This  was  done  to  prevent  surface 
drying  which  could  result  in  detrimental  differential 
stresses  between  a  swollen  wet  core  and  the  dry  outside 
surface  [46].  The  room  temperature  diffusion  coefficient 
for  AS4/3501-6  is  fairly  small,  so  any  additional  water 
absorption  due  to  this  storage  method  should  be  minimal. 

All  specimens  were  usually  tested  within  a  week  of 
conditioning.  Strain  gages  were  applied  after  surface 
drying  and  immediately  prior  to  testing. 

100  grit  sandpaper  tabs  were  bonded  onto  the  gripped 
portions  of  the  specimens  with  a  cyanoacrolate  adhesive. 

A  6.35  mm  thick  specimen  with  sandpaper  tabs  is  shown 
mounted  in  the  compression  fixture  in  Figure  6.  Initial 
specimen  testing  without  the  sandpaper  resulted  in  a 
disproportionate  number  of  end  failures  at  low  strengths. 
The  sandpaper  was  added  to  obtain  a  partial  shear  transfer 
of  the  load  by  increasing  the  frictional  force  between  the 
specimen  and  the  fixture.  Similar  difficulties  with  end 
failures  were  also  encountered  by  Reference  19  who 
concluded  that  20%  of  the  axial  load  had  to  be  transferred 
by  shear  to  avoid  end  failures. 

An  estimate  of  the  frictional  effects  due  to  the 
sandpaper  was  obtained  in  the  following  manner.  The  force 
required  to  push  the  specimen  through  the  fixture  with  and 
without  sandpaper  was  recorded  (see  Figure  13) .  The 
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coefficient  of  friction  without  sandpaper  was  .15  and  the 
coefficient  of  friction  with  sandpaper  was  .08.  This 
means  that  the  per  cent  load  transferred  by  shear  was  48% 
without  sandpaper  and  25%  with  sandpaper.  It  is  apparent 
from  these  measurements  that  the  sandpaper  did  not  increase 
the  frictional  force  between  the  fixture  and  the  specimen. 
However,  the  use  of  sandpaper  eliminated  the  low  strength 
failures  possibly  by  reducing  the  stress  concentration 
where  the  specimen  exited  the  steel  clamping  blocks. 
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ENVIRONMENTAL  CONDITIONING 

The  thick  specimens  for  moisture  conditioning  were 
exposed  to  distilled  water  at  100°  C  (212°  F)  for  29 
days  using  the  experimental  arrangement  illustrated  in 
Figure  14.  The  temperature  of  100°  C  (212°  F)  was  chosen 
so  that  significant  water  absorption  could  be  induced  for 
the  thicker  specimens  in  a  reasonable  amount  of  time. 

Each  specimen  was  surface  dried  and  weighed  after 
conditioning.  SEM  examination  of  a  few  thick  specimens 
following  conditioning  showed  no  evidence  of 
microcracking.  The  average  per  cent  weight  gain  (M)  due 
to  moisture  absorption  for  the  thick  specimens  was  1.86%. 
The  per  cent  weight  gain  following  one  month's  exposure  to 
100°  c  distilled  water  was  predicted  using  Equations  (12) 
and  (13).  Solving  for  M  in  these  equations  required  the 
diffusion  coefficient  (D) ,  the  initial  moisture  content 
(M^ )  and  the  maximum  composite  moisture  content  (Mm) . 

The  diffusion  coefficient  (D)  in  Equation  (13)  was 
determined  as  shown  in  Appendix  B.  The  magnitude  of  this 
coefficient  changes  with  the  inclusion  of  edge  effects 
which  occur  when  the  thickness  of  the  specimen  becomes  a 
significant  fraction  of  its  in-plane  dimensions.  The 
diffusion  of  moisture  through  all  six  sides  of  a 
rectangular  coupon  must  then  be  considered  [38].  Since 
the  cut  edges  of  the  specimens  in  this  study  were  not 
sealed  against  moisture  absorption,  edge  effects  had  to  be 


31 


included  in  this  analysis. 

The  specimens  were  oven  dried  prior  to  conditioning  to 
eliminate  the  initial  moisture  content.  The  maximum 
possible  moisture  content  in  a  composite  laminate  is  the 
product  of  its  resin  weight  fraction  and  the  saturation 
moisture  content  of  the  neat  resin  [41].  Since  the  thick 
specimens  had  a  resin  weight  fraction  of  30%,  and  the 
reported  values  of  the  saturation  moisture  content  of  the 
Hercules  3501-6  neat  resin  range  from  5.5%  to  6.9C 
[41,54],  it  is  evident  that  the  maximum  moisture  content 
of  the  thick  laminate  ranges  between  1.65%  and  2.09%.  The 
measured  moisture  content  of  1.86%  fell  within  this  range 
indicating  that  the  thick  laminates  were  probably 
saturated. 

Ideally,  comparison  of  the  wet  thick  and  thin  laminate 
compression  behavior  should  be  based  on  b: th  laminates 
having  the  same  moisture  content.  However,  the  resin 
weight  in  the  thin  laminates  was  27%  compared  to  the  30% 
resin  weight  in  the  thick  laminates.  Following  the 
process  previously  described,  the  saturation  moisture 
content  of  the  thin  laminates  would  be  less  than  that  of 
the  thick  laminates,  and  range  from  1.5%  to  1.9%. 

Actually,  the  thin  specimens  exhibited  weight  gains  of 
1.3%  after  five  days  of  exposure  to  100°  C  distilled 
water.  These  weight  gains  did  not  change  significantly 
even  with  continued  exposure  for  a  total  of  seventeen 
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days . 

The  measured  saturation  moisture  content  in  the  thin 
laminates  was  slightly  below  that  predicted  using  Equation 
(11).  Several  factors  must  be  considered  when  comparing 
predicted  and  measured  moisture  absorption  values.  First, 
the  predicted  values  require  accurate  knowledge  of  the 
material  diffusion  coefficient.  It  is  not  uncommon  for 
composite  diffusion  coefficients  reported  in  the  literature 
to  be  based  on  data  with  20%  or  more  scatter  [38,55]. 
Secondly,  the  resin  weight  fraction  and  composite  void 
content  must  be  accurately  determined.  Both  of  these 
values  are  typically  determined  by  acid  digestion 
techniques,  which  was  the  method  used  in  this  study. 

The  first  step  in  acid  digestion  is  to  measure  the 
volume  of  water  displaced  by  the  material  sample.  The 
resin  is  then  separated  from  the  fiber  by  soaking  the 
composite  in  nitric  acid.  This  step  allows  the  resin  and 
fiber  weight  fractions  to  be  measured.  If  the  fiber  and 
resin  densities  are  known,  their  respective  volume 
fractions  can  be  calculated.  The  difference  between  the 
sum  of  the  resin  and  fiber  volume  fractions  and  the  total 
volume  of  water  displaced  by  the  sample  is  the  volume  void 
content  in  the  composite  sample.  The  composite  void 
content  obtained  with  this  method  depends  on  the  accuracy 
of  the  fiber  and  resin  densities,  and  on  the  precision  of 
the  measurement  of  the  volume  of  water  displaced  by  the 
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composite  laminate.  Comparison  of  void  contents  measured 
by  acid  digestion  with  those  from  image  processing 
techniques  showed  that  acid  digestion  methods  consistently 
underestimated  the  volume  of  voids  [8].  If  there  is  a 
large  percentage  of  trapped  voids  in  the  resin  in  the 
composite,  more  moisture  could  be  absorbed  than  would  be 
expected  based  on  neat  resin  data  where  the  resin  was 
presumably  void  free. 

All  the  specimens  will  be  considered  saturated  in  this 
study.  This  assumes  that  the  differences  in  moisture 
absorption  between  the  thick  and  thin  specimens  are  due  to 
their  different  resin  contents.  In  addition,  the  thin 
specimens  absorbed  less  moisture  than  predicted,  but  since 
they  showed  no  appreciable  change  in  moisture  content 
between  five  and  seventeen  days  of  moisture  exposure,  it 
is  not  unreasonable  to  assume  these  specimens  were 
saturated. 
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COMPRESSION  TEST  FIXTURE 


An  end  loading  compression  fixture  developed  at  the 
David  Taylor  Research  Center  by  Camponeschi  [9]  was 
selected  for  compression  loading  based  on  its  capacity  for 
testing  6.35  mm  (0.25  in.)  thick  laminates,  and  for  the 
simple  specimen  geometry  required.  The  fixture,  shown  in 
Figure  6,  and  illustrated  schematically  in  Figure  15, 
consists  of  two  hardened  steel  clamping  blocks  each  held 
together  by  four  M6.35  x  1.27  screws.  The  blocks  are  76.2 
mm  x  63.5  mm  x  44.5  mm  (3.0  in.  x  2.5  in.  x  1.75  in.). 
Guidepost  slots  with  linear  bearings  are  employed  to 
maintain  the  fixture  alignment.  Preliminary  trials  with 
the  fixture  showed  that  the  guideposts  were  needed  only  for 
the  thin  specimens.  The  thick  specimens  were  easily 
assembled  in  the  fixture  with  no  guideposts,  and  back-to- 
back  strain  gages  on  these  specimens  during  testing  showed 
no  evidence  of  bending.  Guideposts  were  required  for 
assembly  of  the  fixture  with  the  thin  specimens.  Without 
guideposts,  failures  occurred  during  assembly  as  the  second 
clamping  block  was  mounted  and  tightened  on  the  specimen. 
This  caused  bending  and  failure  of  the  specimen. 
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TEST  PROCEDURE 


Each  specimen  was  mounted  in  the  fixture  with  the  M6.35 
x  1.27  bolts  tightened  to  a  20.25  J  torque.  If  the  bolts 
were  only  hand  tightened,  transverse  brooming  at  the 
specimen  ends  almost  always  occurred  at  relatively  low 
loads.  The  load  was  applied  at  a  rate  of  1.27  mm/mm/min 
(0.05  in. /in. /min. )  using  a  60  kip  Satec  universal 
machine,  with  a  spherically  seated  compression  platen  for 
alignment  purposes.  At  least  ten  specimens  in  each  of  the 
four  series  were  tested  to  failure. 
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RESULTS 


The  mean  dry  thick  compression  strength  was  1055  MPa  + 
19  MPa  (151  ksi  +  3  ksi)  compared  to  1090  MPa  +  70  MPa 
(156  ksi  +  10  ksi)  mean  strength  of  the  dry  thin 
specimens.  The  reported  ranges  for  these  mean  strengths 
were  calculated  using  a  Student  t  distribution  to 
compensate  for  the  small  sample  size.  They  indicate  that 
the  probability  that  the  population  mean  will  be  located 
within  this  range  is  95%  [56].  Statistically,  there  was 
no  significant  difference  between  the  mean  strengths  of 
the  dry  thick  and  thin  specimens.  The  wet  strengths  of 
each  set  of  specimens  were  significantly  lower  than  their 
dry  strengths.  The  mean  thick  wet  strength  was  740  MPa  + 
52  MPa  (106  ksi  +  7  ksi).  The  mean  thin  wet  strength  was 
957  MPa  +  75  MPa  (137  ksi  +  11  ksi) .  The  compression 
strength  data  for  the  thick  and  thin  specimens  is  shown  in 
Tables  3  and  4. 

Stress-strain  curves  for  the  dry  and  wet  specimens  are 
shown  in  Figures  16  through  19.  Due  to  the  small  sample 
size,  no  statistical  analyses  were  performed  on  these 
populations.  The  dry  specimens  (both  thick  and  thin)  had 
initial  modulus  values  of  80  GPa  (11.5  Msi),  as  did  the 
thick  wet  specimens.  The  stress-strain  data  for  the  thin 
wet  specimens,  presented  in  Figure  19,  has  significant 
scatter  with  most  of  the  initial  modulus  values  at  98  GPa 
(14  Msi) .  All  the  observed  stress-strain  curves  showed 
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slight  nonlinearities  typically  beginning  at  40-50%  of  the 
ultimate  compressive  load. 

Typical  thick  and  thin  failure  modes  for  both  the  dry 
and  wet  specimens  are  shown  in  Figures  20  through  23. 
Figures  20  and  21  show  little  difference  in  the  gross 
failure  mode  between  the  thick  dry  and  wet  specimens. 

Each  was  characterized  by  delamination  and  angular 
cracking  across  the  specimen  thickness  to  the  midplane. 
Overall,  delaminations  were  greater  in  the  wet  specimens. 
They  occurred  on  the  outside  surfaces  and  throughout  the 
thickness  of  the  specimen,  often  extending  from  the  gage 
section  to  the  gripped  area  of  the  specimen.  The  crack 
angles  across  the  thickness  ranged  from  27°  to  30°  for 
both  dry  and  wet  specimens,  as  measured  from  an  axis 
perpendicular  to  the  longitudinal  axis  of  the  specimen. 
Occasionally,  biangular  cracking  across  the  specimen 
thickness  was  observed  with  the  mirror  image  cracks 
meeting  in  a  peak  at  the  specimen  midplane.  Failures 
also  occurred  on  some  specimens  at  the  ends  at 
strengths  equivalent  to  those  specimens  that  failed  in 
the  gage  section.  The  failures  at  the  ends  usually 
exhibited  multiple  cracks  through  the  thickness,  again 
at  angles  of  27°  to  30°  for  both  the  wet  and  dry  thick 
laminates . 

The  failed  thin  wet  and  dry  specimens,  pictured  in 
Figures  22  and  23,  also  appeared  very  similar  with 
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slightly  more  delamination  seen  in  the  wet  specimens. 
Direct  comparison  to  their  thick  counterparts  is 
difficult  due  to  the  greater  amount  of  damage  at  failure 
in  the  thin  specimens.  The  thin  specimens  never 
exhibited  clean  angular  cracks  through  the  thickness  as 
observed  in  the  thick  specimens.  Scattered  cracks  across 
the  90°  midplane  plies  were  observed  in  some  of  the  thin 
wet  specimens  as  illustrated  in  Figure  23.  Occasionally, 
these  cracks  were  connected  to  a  delamination  between  the 
90°  and  0°  plies.  Such  cracking  was  not  observed  in  the 
thin  dry  specimens. 
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ANALYSIS 


GLOBAL  FINITE  ELEMENT  ANALYSIS 
MESH  AND  BOUNDARY/LOAD  CONDITIONS 

Three  dimensional  linear  finite  element  analyses  of 
the  dry  and  wet  specimens  were  performed  to  evaluate 
changes  in  the  stress  state  from  swelling  strains  and 
matrix  softening  due  to  moisture  absorption.  The  finite 
element  code  chosen  for  the  analysis  was  ABAQUS  written  by 
Hibbitt,  Karlsson,  and  Sorensen.  It  is  a  general  purpose 
program  featuring  a  number  of  capabilities  that  include 
static  stress,  dynamic  stress,  heat  transfer  and  buckling 
analysis  options.  Non-linearities  may  be  introduced  in 
the  stress  analysis  with  non-linear  constitutive  relations 
that  describe  the  material  (isotropic  materials),  or  by 
geometric  or  boundary  (contact)  non-linearities  [57]. 

The  finite  element  model  of  the  compression  specimen 
and  test  fixture  was  simplified  to  minimize  computational 
costs.  Symmetry  about  the  x-y,  the  y-z  and  the  x-z  planes 
reduced  the  model  to  one-eighth  of  the  actual  specimen 
size  as  shown  in  Figure  24.  The  finite  element  mesh  for 
this  model,  shown  in  Figure  25,  consisted  of  368  20-node 
isoparametric  brick  elements.  As  a  first  approximation, 
clamping  of  the  specimen  by  the  fixture  was  modeled  by 
setting  the  nodal  displacements  in  the  y  direction  equal 
to  zero  in  the  shaded  region  of  Figure  26.  The 
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compression  load  on  the  specimen  was  applied  as 
concentrated  nodal  forces  to  the  nodes  of  surface  A. 

COMPOSITE  MATERIAL  REPRESENTATION 

Ideally,  a  finite  element  model  of  a  composite 
laminate  would  be  constructed  on  a  ply  by  ply  basis.  The 
elements  of  each  ply  could  then  be  represented  as  separate 
orthotropic  materials  with  their  material  properties 
depending  on  the  orientation  of  the  ply  within  the 
laminate.  This  approach  requires  that  each  ply  have  at 
least  one  element  through  its  thickness  which  forces  a 
large  number  of  elements  to  be  used  in  the  other 
directions  to  keep  element  aspect  ratios  reasonable.  For 
example,  a  graphite/epoxy  ply  is  typically  0.127  mm  (0.005 
in.)  thick,  so  that  a  73.00  mm  x  6.35  mm  x  3.18  mm  (2.875 
in.  x  .25  in.  x  .125  in.)  model  would  require  over  700,000 
3-D  elements  if  the  aspect  ratio  of  the  elements  is  1:1. 
Such  an  approach  is  currently  unrealistic  for  thick 
laminates  as  the  large  number  of  plies  would  make  the 
number  of  elements  too  large  for  a  3-D  analysis  to  be 
conducted  at  a  reasonable  cost.  Consequently,  in  this 
analysis,  "smeared"  orthotropic  properties  were  used.  It 
was  assumed  that  the  "smeared"  material  would  have  the 
same  global  response  as  the  material  modeled  on  a  ply  by 
ply  basis  [ 58 ] . 
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MOISTURE  EFFECTS 


Two  changes  were  made  in  the  model  to  accommodate  the 
moisture  conditioned  specimens.  First,  the  smeared 
material  properties  were  adjusted  to  account  for  the  loss 
of  stiffness  in  the  unidirectional  transverse  modulus  due 
to  moisture.  This  adjustment  involved  a  reduction  in  the 
transverse  Young's  modulus  (E2  and  E3)  of  dry  AS4/3501-6 
[59]  by  12%  to  get  the  unidirectional  wet  lamina 
properties.  The  reduction  of  12%  was  selected  based  on 
experimental  data  reported  by  Chi-Hung  and  Springer  [60] 
for  AS4/3501.  The  (02/90)s  cross  ply  laminate  smeared 
material  properties  were  then  obtained  using  STIF3D  which 
is  PC-based  program  developed  at  the  University  of 
Delaware  [61].  The  effective  homogeneous  constitutive 
equations  used  in  STIF3D  were  originally  reported  by 
Pagano  [62],  and  include  the  effects  of  bending- 
extensional  coupling.  The  resultant  3-D  unidirectional 
and  cross  ply  material  properties  for  the  dry  and  wet 
material  are  shown  in  Table  5.  Since  the  cross  ply 
laminate  is  dominated  by  the  0°  plies,  stiffness 
reductions  in  the  unidirectional  transverse  modulus  did 
not  significantly  affect  the  smeared  material  properties. 

The  second  change  made  in  the  model  to  represent  the 
effects  of  moisture  was  the  incorporation  of  swelling 
strains.  For  a  unidirectional  AS4/3501  lamina,  the 
coefficient  of  moisture  expansion,  p,  (analagous  to  the 
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coefficient  of  thermal  expansion)  as  reported  by  Tsai  [63] 
is  negligible  in  the  fiber  direction.  In  the  transverse 
and  through  thickness  directions  (5  is  0.4  4  %M  The 

moisture  expansion  coefficients  for  a  (0/90)  laminate, 
based  on  classical  laminated  plate  theory  are  then  0.147 
%M  1  in  the  0°  direction,  0.293  %M  ^  in  the  90°  direction 
and  0.44  %M_1  through  the  thickness. 

Multiplication  of  the  coefficient  of  moisture 
expansion  by  the  moisture  concentration  yields  the 
swelling  strain  which  can  then  be  superimposed  on  the 
mechanical  strains  due  to  the  compressive  load. 

Generation  of  the  moisture  concentration  profile  through 
the  specimen  thickness  was  based  on  Fickian  diffusion  and 
corrected  for  edge  effects  using  the  methodology  presented 
by  Springer  [38].  All  calculations  are  shown  in  Appendix 
B.  The  moisture  concentrations  as  a  function  of  position 
through  the  thickness  of  both  the  thick  and  thin  specimens 
are  shown  in  Figures  27  and  28  for  two  exposure  time 
periods.  These  results  show  that  the  distribution  of  the 
moisture  concentration  was  nearly  constant  for  the  thick 
and  thin  specimens. 

Thermal  strains  associated  with  a  temperature 
differential  are  easily  calculated  in  ABAQUS.  Since 
swelling  strains  are  analogous  to  thermal  strains,  the 
moisture  concentration  data  from  Figures  27  and  28  was 
simulated  in  ABAQUS  as  a  temperature  distribution,  and  the 
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moisture  expansion  coefficients  were  simulated  as  thermal 
expansion  coefficients.  Since  the  moisture  gradients  were 
very  small,  and  the  changes  in  the  elastic  constants  were 
moderate,  large  differences  between  the  wet  and  dry 
analyses  were  not  expected. 

RESULTS 

The  distribution  of  the  stresses  (three  normal  and 
three  shear)  showed  only  minor  differences  between  the 
stress  state  existing  in  the  wet  and  dry  specimens.  The 
longitudinal  stress  (o2)  was  typically  two  to  three 

orders  of  magnitude  greater  than  the  remaining  stresses 
(oy»ox,Txy,T12,Tyz)  in  both  the  wet  and  dry  specimens. 

Slight  increases  in  ox  and  oy  occurred  due  to  the 

addition  of  moisture,  but  the  overall  magnitude  of  the 
stresses,  ox  and  oy,  remained  small  (less  than  7  MPa)  . 
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MICROMECHANICAL  ANALYSIS 
BACKGROUND 

In  addition  to  the  global  stress  analysis  conducted 
using  ABAQUS,  a  micromechanic  analysis  was  conducted  to 
compare  the  compression  strength  of  wet  and  dry  specimens. 
As  discussed  in  the  literature  review,  the  compression 
failure  theory  of  Reference  [6]  was  chosen  for  this 
comparison.  Use  of  this  failure  theory  requires  the  shear 
stresses  in  the  matrix  as  a  function  of  the  compressive 
load  on  the  composite  laminate.  The  approach  presented  in 
Reference  [6]  to  determine  the  matrix  shear  stresses  served 
as  a  guideline  for  the  analysis  in  this  study,  and  will- 
therefore,  be  briefly  discussed. 

In  [6],  a  composite  laminate  was  modeled  as  a 
multilayered  Timoshenko  beam  that  permitted  transverse 
deflections  of  each  layer  and  allowed  the  cross  sections 
to  rotate  with  respect  to  the  neutral  axis.  Using  this 
model,  two  cases  were  considered.  The  first  case 
represented  the  entire  laminate  as  alternating  fiber  and 
matrix  layers.  Since  this  was  not  practical  for  thicker 
laminates,  a  second  case  was  evaluated  which  consisted  of 
a  fiber  layer  surrounded  by  two  matrix  half  layers, 
representing  a  repeating  structure  within  the  laminate. 
Expressions  for  the  displacement  of  each  layer  were 
obtained  using  multilayered  Timoshenko  beam  theory.  These 
diplacements  represented  those  due  to  bending  and  shearing 
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of  the  cross-section.  The  axial  displacements  due  to  the 
axial  compressive  load  on  the  beam  were  ignored. 
Displacements  were  converted  to  strains  using  the  strain- 
displacement  equations,  and  Hooke's  law  was  used  to 
convert  the  strains  to  stresses.  A  set  of  linear  second 
order  partial  differential  equations  for  this  problem  were 
then  developed  using  the  equilibrium  relations.  These 
equations  were  expressed  in  finite  difference  form  and 
solved  numerically.  The  nonlinearity  of  the  matrix 
material  was  incorporated  in  the  analysis  by  dividing  the 
matrix  shear  stress-strain  curve  into  a  number  of  linear 
segments.  Loads  were  then  applied'incrementally  with 
variations  in  the  matrix  shear  modulus  occuring  within  the 
load  increments.  An  initial  imperfection  in  the  beam  was 
allowed  and  mathematically  described  by 

wQ=a0{l+sin[  (^f)  (x-J)  ]  }  (23) 

wo=  amplitude  of  imperfection  as  a  function  of 
length  along  the  beam 
a0=  half  the  maximum  amplitude  of  the 
imperfection 

L=  length  of  the  beam  equivalent  to  half 
the  wavelength  of  the  imperfection 

The  solution  of  the  governing  set  of  differential 
equations  yielded  displacements  which  could  be  used  to 
calculate  the  shear  stresses  in  each  layer. 

The  shear  stresses  obtained  with  the  preceding 
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solution  depended  on  which  case  (entire  laminate  or 
repeating  element)  was  modeled.  The  repeating  element 
model  typically  predicted  shear  stresses  that  were  15-20% 
lower  than  those  predicted  with  a  model  of  the  entire 
laminate  [33].  The  difference  in  the  predicted  shear 
stresses  between  the  two  models  increased  with  a  decrease 
in  the  fiber  to  matrix  longitudinal  modulus  ratio. 

Two  failure  theories  were  developed  based  on  these 
calculated  shear  stresses.  If  the  shear  stress  in  the 
matrix  layer  exceeded  the  interlaminar  shear  strength  of 
the  material,  failure  would  occur  by  delamination.  The 
second  mode  of  failure  was  shear  buckling.  Assuming  that 
the  lamina  length-to-width  ratio  was  small  enough  to 
preclude  Euler  buckling,  Davis  postulated  that  the  change 
in  strain  energy  from  the  initial  to  the  final  position 
was  due  solely  to  shearing  strains.  Equation  (5)  is  the 
resultant  expression  for  the  critical  axial  compressive 
load  leading  to  shear  instability. 

The  compression  strength  of  the  composite  is  based  on 
the  shear  moduli  of  the  fiber  and  matrix  and  their 
respective  volume  fractions.  Since  the  fiber  shear 
modulus  and  volume  fraction  do  not  change  with  moisture 
absorption,  the  matrix  shear  modulus  will  be  the  critical 
parameter.  Accurate  use  of  this  failure  theory  is  then 
dependent  on  two  factors.  First,  the  matrix  shear  stress- 
strain  curve  must  be  known  as  a  function  of  moisture.  And 
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secondly,  the  correct  matrix  shear  modulus  corresponding 
to  the  shear  stress-strain  state  of  the  matrix  must  be 
chosen. 

MESH  AND  BOUNDARY  CONDITIONS 

In  this  analysis,  the  shear  stresses  in  the  matrix  of 
a  composite  loaded  in  compression  were  determined  with 
finite  element  techniques  (ABAQUS  program)  versus  finite 
difference  methods.  Axial  displacements  due  to  axial 
compression  were  included  in  the  analysis.  Two  models  of 
the  composite  laminate  were  considered.  In  the  first 
case,  a  repeating  element  within  the  laminate  was  modeled 
as  a  fiber  plate  surrounded  by  two  half  layers  of  matrix. 
The  fiber  and  matrix  layer  thicknesses  were  proportional 
to  their  respective  volume  fractions  within  the  ply.  The 
total  thickness  of  the  beam  was  equal  to  one  ply  thickness 
0.127  mm  (0.0055  in.).  This  model  was  very  similar  to 
Davis'  repeating  element  case  for  modeling  boron/epoxy 
composites.  However,  the  diameter  of  a  boron  fiber  is 
approximately  fourteen  times  that  of  a  graphite  fiber,  so 
modeling  the  fibers  in  a  boron/ epoxy  ply  as  a  single  fiber 
plate  is  reasonable.  Since  the  graphite  fibers  are  so 
small  in  diameter,  many  fibers  will  be  stacked  through  the 
ply  thickness,  so  the  single  fiber  plate  model  might  not 
reflect  interactive  effects  between  fibers.  Therefore,  a 
second  model  was  also  included  in  this  analysis  which 
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consisted  of  two  fiber  plates  separated  by  a  matrix  layer 
and  surrounded  by  matrix  half  layers.  The  total  fiber  and 
matrix  layer  thicknesses  were  again  proportional  to  their 
respective  volume  fractions  within  the  ply.  A  4400 
element  mesh  was  used  for  each  of  these  models  and  is 
shown  in  Figure  29.  The  elements  were  second  order 
(8  node)  plane  strain,  isoparametric  quadrilaterals. 

The  initial  imperfection  in  the  fiber-matrix  beam  was 
generated  using  Equation  (21)  with  two  amplitude-to-length 
ratios,  0.00375  and  0.009.  These  ratios  were  chosen 
somewhat  arbitrarily  as  data  on  initial  imperfections  in 
flat  graphite/epoxy  laminates  is  not  extensively 
available.  The  lower  ratio  represents  amplitudes  measured 
by  Davis  for  boron/epoxy  laminates,  and  Camponeschi  for 
graphite/ epoxy  laminates  [33,64).  Significantly  higher 
ratios  (0.025)  were  used  by  Hyer  and  Telegadas  to 
represent  waves  in  graphite/ epoxy  filament  wound  cylinders 
[7,65].  Imperfections  this  large  would  be  unlikely  in  a 
flat  laminate.  The  ratio  of  0.009  was  chosen  as  a  midway 
point  between  wave  amplitudes  seen  in  flat  plates  and 
filament  wound  cylinders,  and  to  illustrate  the  effects  of 
an  increasing  wave  amplitude  on  the  predicted  shear 
stresses  in  the  matrix. 

The  beam  length  (or  half  wavelength)  was  0.76  mm 
(0.03  in.),  which  is  equivalent  to  half  the  critical 
length  for  Euler  buckling.  This  resulted  in  a  beam 
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length-to-thickness  ratio  of  5.5:1.  No  attempt  was  made  to 
evaluate  the  effects  of  variations  in  wavelength  on  the 
magnitude  of  the  shear  stresses  at  the  fiber-matrix 
interface. 

Uniform  axial  displacements  corresponding  to  1.5% 
strain  for  the  dry  specimens  and  1.2%  strain  for  the  wet 
specimens  were  imposed  on  one  end  of  the  beam  while  the 
other  end  was  held  fixed.  The  model  had  been  sized  to 
avoid  Euler  buckling,  but  the  analysis  was  geometrically 
non-linear  allowing  for  large  displacements  and  rotations 
should  they  occur  due  to  material  property  changes. 
Therefore,  no  symmetry  boundary  conditions  were  used. 

The  top  and  bottom  surfaces  of  the  model  were  free  to 
move  in  the  x  and  y  directions.  In  actuality,  some  type 
of  normal  and  shear  forces  are  exerted  on  these  surfaces 
by  the  plies  above  and  below  the  repeating  element.  The 
boundary  conditions  on  these  surfaces  are  difficult  to 
describe,  particularly  in  the  presence  of  a  wavy  ply. 

Davis  [33]  fixed  the  axial  displacement  of  the  nodes  on 
these  surfaces.  Hyer  and  Telegadas  [65]  modeled  several 
layers  above  and  below  a  wavy  ply,  and  then  used  an 
elasticity  solution  for  boundary  conditions  away  from  the 
wave. 

Neither  of  these  conditions  were  appropriate  for  the 
model  in  the  current  analysis  which  differed  slightly  from 
the  preceding  models.  Preliminary  trials  with  the  current 
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model  showed  a  reduction  in  shear  stress  with  the  addition 
of  any  normal  forces  on  the  top  and  bottom  surfaces  of  the 
model.  Therefore,  no  constraints  were  placed  on  the  nodes 
on  these  surfaces  in  an  attempt  to  simulate  a  worst  case 
scenario. 


MATERIAL  MODEL 

The  fiber  was  modeled  as  a  transversely  isotropic 
material  with  the  properties  shown  in  Table  1.  A  three 
parameter  Richard-Blacklock  equation  was  used  to  describe 
the  nonlinear  shear  behavior  of  the  matrix  [48].  The 
relationship  established  between  shear  stress  and  shear 
strain  is  shown  in  Equation  (22). 


(24) 


G=  initial  shear  modulus 

n=  curvature  parameter 

to=  asymptotic  shear  stress  value 


Adams  developed  expressions  for  the  parameters  in  Equation 
(24)  as  functions  of  temperature  and  moisture  [49].  The 
matrix  shear  stress-strain  curves  generated  using  this 
model  are  shown  in  Figure  30  for  the  dry  and  wet  (matrix 
moisture  contents  corresponding  to  1.3%  and  1.9%  moisture 
contents  by  weight  in  the  composite  laminate)  material. 

If  the  shear  stress  in  the  matrix  is  known  as  a  function 
of  axial  compressive  load,  the  proper  modulus  can  be 
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chosen  from  these  curves  for  use  in  Equation  (5) .  The 
finite  element  analyses  were  performed  with  the  matrix 
material  behaving  linear  in  shear.  An  attempt  was  made  to 
bound  the  solution  by  using  three  values  for  the  matrix 
shear  modulus  (initial  dry,  intermediate,  and  final  wet) 
for  each  amplitude  to  wavelength  ratio.  These  modulus 
values  were  1.6  GPa  (.23E6  psi) ,  0.99  GPa  (.1415E6  psi)  , 
and  0.37  GPa  (0.53E5  psi),  respectively.  This  covered  the 
range  of  matrix  shear  modulus  values  from  its  stiffest 
value  (initial  dry  modulus)  to  its  softest  value  (final 
wet  modulus) .  The  analysis  with  a  uniform  axial 
displacement  corresponding  to  a  strain  of  T.2%  (wet 
specimens)  was  performed  using  only  the  intermediate 
matrix  modulus  value  for  both  amplitude  to  length  ratios. 

RESULTS 

Stresses  in  the  global  coordinate  system  were  obtained 
at  the  integration  points  from  the  finite  element 
analysis.  The  through  thickness  shear  stresses  were 
expected  to  be  continuous  based  on  classical  laminated 
plate  theory  which  assumes  equal  and  opposite  traction 
boundary  conditions  between  plies  for  the  out-of-plane 
stresses  [66].  However,  this  required  a  rotation  of  the 
stresses  from  the  global  coordinate  system  to  a  local 
coordinate  system  normal  and  tangential  to  the  element 
edges  [7,65].  The  highest  matrix  shear  stresses  were 
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expected  to  occur  at  the  fiber-matrix  interface. 

The  maximum  interfacial  shear  stresses  occurred  at 
points  one-fourth  and  three-fourths  the  length  of  the  beam 
for  both  the  single  fiber  and  double  fiber  plate  models. 
The  matrix  and  fiber  shear  stresses  at  these  quarter 
points  A  (Figure  29  a)  and  b) )  in  the  single  and  double 
fiber  plate  models  were  very  similar.  However,  in  the 
double  fiber  plate  model,  the  matrix  shear  stresses  at 
point  B  in  Figure  29  b)  were  typically  8  to  22  times 
greater  than  the  matrix  shear  stresses  at  point  A.  These 
results  indicated  that  interactive  effects  resulting  from 
many  fibers  stacked  through  the  ply  thickness  could  have 
a  significant  effect  on  the  predicted  matrix  shear 
stresses.  Therefore,  the  double  fiber  plate  model  was 
chosen  as  the  prototype  for  this  study.  The  maximum 
matrix  shear  stresses  in  this  model  occurred  along  the  C-D 
interface  shown  in  Figure  29. 

The  maximum  shear  stresses  (at  the  integration  points) 
in  the  fiber  and  matrix  elements  on  either  side  of  the  C-D 
interface  are  shown  in  Figures  31  and  32  for  both 
amplitude-to-wavelength  ratios.  To  be  conservative,  an 
average  of  the  shear  stresses  from  the  fiber  and  matrix 
integration  points  on  either  side  of  the  interface  was 
chosen  to  represent  the  maximum  shear  stresses  in  the 
matrix  at  the  interface.  Additional  mesh  refinement  at 
the  interface  would  be  required  to  predict  the  interfacial 
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shear  stresses  more  exactly.  Table  6  lists  the  maximum 
shear  stresses  from  each  of  the  analyses  which  will  be 
used  to  select  tangent  shear  moduli  from  the  matrix  shear 
stress-strain  curves  in  Figure  30.  These  tangent  shear 
moduli  will  then  be  used  in  Equation  (5)  to  predict  the 
compression  strength  of  the  graphite/ epoxy  laminates. 
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DISCUSSION  AND  CONCLUSIONS 


Experiment 

End  loading  compression  techniques  have  given  lower 
compression  strengths  for  composites  as  compared  to 
shear  loading  methods  (IITRI  and  Celanese  compression 
tests,  for  example) .  Since  the  laminate  strengths  in 
this  study  were  generated  using  end  loading  techniques, 
it  was  important  to  establish  some  confidence  in  the 
test  method  by  comparing  these  results  with  the  strength 
determined  by  other  methods. 

The  mean  laminate  compressive  strengths  were  1055  MPa 
(151  ksi)  and  1090  MPa  (156  ksi)  for  the  dry  thick  and 
thin  laminates,  respectively.  Using  classical  laminated 
plate  theory  [62],  a  unidirectional  longitudinal 
compressive  modulus  of  119  GPa  (17E6  psi)  [9],  and  a 
unidirectional  transverse  compressive  modulus  of  10  GPa 
(1.5E6  psi)  [66],  the  critical  stress  was  1484  MPa  (212.5 
ksi)  and  occurred  in  the  0°  plies  of  the  (02/90)s 
laminate.  Compressive  unidirectional  lamina  strengths  for 
AS4/3501-6  generated  using  the  IITRI  compression  fixture 
range  from  1397  to  1537  MPa  (200  to  220  ksi) 

[21,23,28,50].  Since  the  stresses  developed  in  the  0° 
plies  in  this  study  were  within  this  range,  it  was 
concluded  that  the  test  method  developed  was  satisfactory. 

The  literature  review  indicated  that  a  drop  (20%  or 
less)  in  compression  strength  was  expected  after  moisture 
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conditioning  to  saturation.  We  found  that  the  wet 
laminate  strengths  were  11%  and  30%  lower  than  the  dry 
laminate  strengths  for  the  thin  and  thick  specimens, 
respectively.  The  larger  loss  of  strength  in  the  thick 
laminates  may  be  due  to  the  greater  amount  of  moisture 
absorbed  by  these  specimens.  Linear  extrapolation  of  the 
strength  of  a  specimen  with  1.3%  moisture  to  one  with 
1.86%  moisture  (observed  in  the  thick  laminates)  results 
in  a  drop  in  the  compresssion  strength  to  678  MPa  (97 
ksi) .  This  prediction  is  consistent  with  the  actual 
strengths  observed  in  the  wet  thick  laminates. 

Moisture  had  little  effect  on  the  fiber  dominated 
longitudinal  modulus  between  0%  and  0.2%  strain  in  the 
thick  laminates.  The  thin  laminates  showed  an  increase  in 
longitudinal  modulus  after  moisture  conditioning  which  is 
not  consistent  with  any  cases  reported  in  the  literature. 

It  appears  that  the  outer  plies  of  the  thin  wet  laminate 
buckled,  and  the  strain  gage  sensed  tensile  strains  due  to 
bending.  These  tensile  strains,  in  the  presence  of  the 
compressive  strain  due  to  the  longitudinal  load,  lead  to  a 
smaller  combined  strain  measurement  (and  higher  modulus) 
as  compared  to  the  case  with  no  outer  ply  buckling. 

A  comparison  in  the  gross  failure  mode  of  the  wet  and 
dry  specimens  in  Figures  20  through  23  for  the  thick  and 
thin  laminates  showed  little  difference  except  for  slightly 
greater  delaminations  in  the  wet  specimens.  These  results 
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are  similar  to  those  reported  by  Grimes  and  Adams  [50]. 

The  catastrophic  nature  of  the  final  failure  destroyed  any 
clues  as  to  the  failure  mechanisms  immediately  preceding 
final  failure  that  might  be  determined  by  microscopic 
examination.  It  is  possible  that  different  mechanisms 
initiated  failure  in  the  wet  and  dry  laminates,  but  that 
the  gross  features  of  the  failure  appeared  the  same. 

The  thick  specimens  typically  exhibited  a  combination 
of  delamination  and  angular  cracking  through  the  laminate 
thickness.  Similar  failure  modes  under  compressive 
loading  have  been  observed  in  dry  graphite/epoxy  by  Smoot 
[23],  Garala  [25]  and  Camponeschi  [64],  The  crack  angles 
through  the  thickness  ranged  from  27°  to  30°  for  both  the 
wet  and  dry  thick  laminates.  These  angles  are 
consistent  with  the  26°  kink  band  angles  reported  by 
Hahn  [1]  for  graphite/ epoxy  composites. 

The  angularly  cracked  surfaces  were  not  observed 
following  failure  in  the  thin  specimens.  The  majority 
of  these  specimens  had  significant  damage  in  the 
failure  zone  making  identification  of  the  failure  mode 
difficult.  The  wet  specimens  exhibited  more 
delamination,  and  a  few  showed  significant  cracking 
across  the  90°  midplane  plies.  This  may  be  a  result 
of  reduced  matrix  properties  due  to  moisture 
absorption. 
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Theory 

The  differences  in  stress  states  predicted  by  the  3-D 
finite  element  analysis  for  the  wet  and  dry  laminates  did 
not  explain  the  large  differences  in  strength  measured  in 
the  laboratory.  This  result  was  expected  because  the 
moisture  gradients  modeled  in  the  analysis  were  small,  and 
large  differences  in  the  stresses  between  the  swollen  wet 
layers  and  the  dry  layers  did  not  develop. 

Since  moisture  absorption  in  a  composite  material 
occurs  in  the  matrix,  it  is  not  unreasonable  to  assume 
that  changes  in  matrix  behavior  with  moisture  are  linked 
to  the  compression  strength  differences  of  the  dry  and  wet 
laminates.  A  finite  element  analysis  on  the  constituent 
level  was  performed  t~»  determine  the  shear  stresses  at  the 
fiber-matrix  interface  in  a  composite  ply  loaded  in 
compression.  Initial  waviness  of  the  ply  was  incorporated 
in  the  model.  The  maximum  interfacial  shear  stress  from 
the  analysis  was  used  to  select  a  tangent  modulus  from  the 
wet  and  dry  matrix  shear  stress-strain  curves.  The 
appropriate  tangent  shear  modulus  was  then  used  in  Davis' 
[6]  compression  failure  theory  (Equation  (5))  based  on 
shear  instability. 

The  micromechanical  analysis  showed  maximum  matrix 
shear  stresses  ranging  from  34  MPa  (4.8  ksi)  to  65  MPa 
(9.3  ksi)  for  an  initial  imperfection  amplitude  to 
wave.lengt:  ratio  of  0.00375,  and  from  70  MPa  (10.0 


58 


ksi)  to  140  MPa  (20.0  ksi)  for  a  ratio  of  0.009  (See 
Table  6) .  The  upper  bound  of  each  range  occurred  when 
the  initial  dry  matrix  shear  modulus  was  used  in  the 
analysis,  and  the  lower  bound  occurred  with  the  use  of 
the  final  wet  modulus.  Since  the  matrix  shear  stress- 
strain  behavior  is  highly  nonlinear,  use  of  these 
moduli  might  overestimate  or  underestimate  the  actual 
matrix  shear  stresses.  Therefore,  analyses  were  also 
run  with  an  intermediate  matrix  shear  modulus. 

Predicted  matrix  shear  stresses  then  ranged  from  54 
MPa  (7.8  ksi)  to  112  MPa  (16.0  ksi)  for  the  1.5% 
strain  case  and  from  44  MPa  (6.2  ksi)  to  91  MPa  (13.0 
ksi)  for  the  1.2%  case. 

The  upper  bound  stress  from  each  of  these  ranges 
exceeds  the  reported  maximum  interlaminar  shear  strength 
of  91  MPa  (13.0  ksi)  for  the  dry  material  and  70  MPa 
(10.0  ksi)  [67]  for  the  wet  material.  A  wave  amplitude 
to  length  ratio  of  0.009  was  used  to  predict  the  upper 
bound  stresses.  This  ratio  may  be  too  large  for 
autoclave  cured  laminates.  The  upper  bound  stress  of 
each  range  was  reduced  to  the  appropriate  (dry  or  wet) 
interlaminar  shear  strength.  However,  predicted  matrix 
shear  stresses  were  of  sufficient  magnitude  to  indicate 
that  delamination  was  also  a  possible  failure  mode  for 
these  laminates.  Experimental  observations  of  combined 
failure  modes  (delamination  and  shear  cracking  through 
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the  specimen  thickness)  support  this  conclusion. 

The  failure  stresses  determined  experimentally  are 
compared  with  the  stresses  predicted  from  Equation  (5)  in 
Table  7.  Two  predicted  stresses  are  shown.  The  first 
uses  a  matrix  tangent  shear  modulus  corresponding  to  the 
lower  bound  matrix  shear  stresses,  and  the  second  uses  a 
matrix  shear  modulus  at  the  interlaminar  shear  strength  of 
the  material.  The  second  prediction  correlates  best  with 
the  experimental  results. 

CONCLUSIONS 

The  compression  behavior  of  thick  (6.35  mm)  (02/90)s 
AS4/3501-6  laminates  was  compared  to  that  of  thin  (1.6  mm) 
section  laminates  in  both  dry  and  moisture  saturated 
states.  Experimental  results  showed  no  statistical 
difference  in  strength  of  the  dry  material  of  either 
thickness.  However,  there  was  a  statistically  significant 
difference  in  the  strengths  of  moisture  conditioned 
specimens  of  different  thicknesses.  The  thick  specimens 
absorbed  significantly  more  water  than  the  thin  specimens, 
and  consequently,  showed  a  greater  loss  in  compression 
strength.  This  experimental  observation  was  substantiated 
with  a  stress  analysis  that  accounted  for  the  reduced 
matrix  properties  of  the  thicker  specimens  associated  with 
the  higher  moisture  content. 

Slightly  more  delamination  was  observed  in  the  wet 
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specimens  of  both  thicknesses.  The  thick  laminates 
exhibited  a  combination  of  delamination  and  angled 
cracks  of  27°  to  30°  partially  across  the  thickness  of 
the  specimen.  Comparison  of  the  failure  modes  of  the 
thick  and  thin  laminates  was  difficult  due  to  the 
extensive  delamination  and  damage  that  occurred  in  the 
thin  laminates.  Well  defined  angular  cracks,  as  seen 
in  the  thick  specimens,  were  not  observed. 

No  change  in  longitudinal  modulus  was  observed  with 
moisture  concentration  or  with  thickness. 

Strength  differences  between  the  dry  and  wet  laminates 
are  not  due  to  moisture  induced  stresses.  It  appears  that 
the  detrimental  influence  of  the  moisture  is  the 
degradation  of  the  matrix  shear  modulus  at  high 
compressive  stresses  with  high  moisture  concentrations. 
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APPENDIX  A 


BUCKLING  LOAD 


The  equation  used  to  calculate  the  global  buckling 
stresses  for  the  compression  specimens  was  derived  by 
Timoshenko  and  included  transverse  shear  effects  [68]. 
Clamped-clamped  end  conditions  were  assumed. 


°cr= 


Ex7H 


12(|)Wl.2gi 


(A-l) 


ocr=  critical  buckling  stress 
Ex=  longitudinal  modulus 
/=  unsupported  length 
t=  specimen  thickness 
Gi2=  transverse  shear  modulus  - 

For  AS4/3501— 6  in  a  (0/90)  layup,  the  longitudinal 

/  s 

modulus  is  84  GPa  (12E6  psi) ,  and  the  transverse  shear 
modulus  is  estimated  at  3 . 5E6  GPa  ( . 5E6  psi).  These 
modulus  values  were  not  changed  for  the  conditioned 
specimens  as  a  large  effect  on  the  longitudinal  modulus 
was  not  expected  for  this  particular  layup.  Safety 
factors  were  based  on  an  estimated  failure  strength  of 
1048  MPa  (150  ksi) . 

Thick  Specimens  3  ccr=  2109  MPa  (302  ksi) 

S.F.=2 . 0 


Thin  Specimens  6  ocr=  1152  MPa  (165  ksi) 

S . F. =1 . 1 
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APPENDIX  B  MOISTURE  CONCENTRATION 


Prediction  of  the  distribution  of  the  moisture 
concentration  through  the  thickness  of  the  specimen 
requires  a  diffusion  coefficient  which  takes  into  account 
edge  effects  and  laminate  orientation.  Accounting  for 
edge  effects  was  necessary  in  this  program  because  the 
specimen  edges  were  not  coated  during  immersion  in  100°  c 
(212°  F)  distilled  water. 

Calculation  of  the  resin  diffusion  coefficient  is  the 
first  step  in  estimating  the  diffusion  coefficient  for  the 
composite  laminate.  An  empirical  expression  for  the 
diffusion  coefficient  of  a  unidirectional  AS4/3501 
composite  with  a  fiber  volume  fraction  of  65%  immersed  in 
distilled  water  is  given  in  Equation  B-l  [38]. 

-7218 

Dc=768e  T  (B-l) 

Dc=  composite  diffusivity  in 

T=  temperature  in  Kelvin 

The  resin  diffusivity  is  given  by 


Dr=  resin  diffusivity  in 
Vf=  fiber  volume  fraction 

For  100°  C  exposure, 
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Dc=  3 . 0286E-6  ^gr  (1.6900E-5  |  -) 

Dr=  3 . 3550E-5  ^  (1.8721E-4 

Springer  derived  the  following  expression  for  the 
diffusion  coefficient  of  a  laminated  coupon  including  edge 
effects  [38]. 

D=  Dfe(l-2>|5)[1+^ 


(B-3) 

p^=  90°-v.  (See  Figure  10) 

j=  1  to  N,  N=  number  of  plies  in  laminate 

For  a  (0/90)  layup,  Equation  B-3  reduces  to: 

2  S 

Thick  Specimens 

Zhj=  6.35  mm  (0.25  in.) 

Zhjcos^  | (6. 35) 

IhjSin2^  ^(6.35) 

D=.  11DR  ( 1+1 . 624y+  1 . 348^)  2  V,=  .62  (B-4) 

Thin  Specimens 

Zhj=  1.6  mm  (0.062  5  in.) 

ZhJcos20.=  |  (1.6) 

ZhjSin2^  ^(1.6) 

D=  .  07  Dp  ( 1+1 . 839y+  1.480^)2  Vf=.68  (B-5) 
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.  .  h  6.35  h  6.35 

For  the  thick  specimens,  y=  i4  6~  n=i2  .~7 

Dmici=  1  •  123E-5  ^(6.265E-5^?) 

„  ,  ,  ,  ,  .  h  1.6  h  1.6 

For  the  thin  specimens,  j=x36T5  n=i2 7l 

Dmln=  3.418E-6  3^(1.907E-5 

A  computer  program  was  written  to  calculate  C(z,t) 
from  Equation  (9)  in  the  main  text  using  these  diffusion 
coefficients. 
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Table  1  Elastic  Constants 
for 

AS4  Fiber 


Elastic  Constants 

AS4  Fiber 

GPa  (psi) 

El 

237  (34E6) 

E 2 

14  (2E6) 

E3 

1 4  (2E6) 

v  12 

20 

vl3 

.30 

v23 

.30 

G12 

1 4  (2E6) 

G13 

14  (2E6) 

G23 

5  (.77E6) 

Fiber  Radius  7E-6  m  (2.76E-4  in.) 


Table .?  Elastic  Constants 
for 

3501-6  Resin 


Elastic  Constants 

3501-6  Resin 

GPa  (ds\) 

E 

4  (.62E6) 

V 

.34 

G 

1.6  (.23E6) 
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Table  3  Compression  Strength 
for  Thick  (6.35  mm) 
AS4/3501-6  (Cy90)s 
Composites 


S  x  Standard  Deviation 

S7  Standard  Deviation  of  the  Mean 

t(a)  Student's  t  Statistic 
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Table  4  Compression  Strength 
for  Thin  (1.6  mm) 
AS4/3501-6  (02/90). 
Composites 


Standard  Deviation 

O— - 

*  Standard  Deviation  of  the  Mean 
t(a)  Student's  t  Statistic 
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Table  5  Elastic  Constants  for  Dry  and  Wet 
(Saturated,  Room  Temperature) 
AS4/3501 -6  Laminates 
V  =  60% 


Unidirectional 

Dry  Wet 

GPa  (psi)  GPa  (psi) 

(0j/90) 

Dry 

GPa  (psi) 

Laminate 

Wet 

GPa  (psi) 

Ell 

117(16.7E6) 

117(16.7E6) 

82(1 1 .7E6) 

81(11.6E6) 

E22 

10(1.5E6) 

95(1. 3E6) 

46(6.6E6) 

45(6.5E6) 

E33 

10(1 .5E6) 

9.2(1.3E6) 

12(1 .8E6) 

10(1 .5E6) 

vl2 

.33 

.33 

.075 

.067 

vl3 

.33 

.451 

.455 

v23 

.47 

.47 

.47 

.47 

G12 

6(.87E6) 

6(.87E6) 

6(.87E6) 

6(.87E6) 

G13 

6(.87E6) 

6(.87E6) 

5(.73E6) 

5(.73E6) 

G23 

3.8(.55E6) 

3.8(.55E6) 

4.4(.63E6) 

4.4(.63E6) 

1  Fiber  Direction  1  o' Direction 

2  Transverse  to  Fiber  2  90'Direction 

3  Through  Thickness  3  Trhough  Thickness 
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Table  6  Maximum  Matrix  Shear  Stresses 
Due  to  Axial  Compressive  Loads 
In  An  Initailly  Wavy  AS4/3501-6  Ply 


a/L 

g; 

GPa(Msi) 

0.00375 

1.60(0.23) 

0.99(0.14) 

0.37(0.05) 

0.009 

1.60(0.23) 

0.99(0.14) 

0.37(0.05) 

Maximum  Matrix  Shear  Stress 
MPa  (ksi) 


1 ,2%Strain  1 .5%  Strain 


65.0(9.3) 

54.0(7.8) 

34.0(4.8) 

140.0(20.0) 


68.0(9.7) 
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Table  7  Comparison  of  Predicted  and  Measured 
Compression  Strengths 


vF 

vM 

Gp 

GPa(Msi) 

^MT 

MPa(Msi) 

o 

uni 

MPa(ksi) 

Predicted 

a 

0j/90 

MPafksil 

Measured 

o 

02/90 

MPafksi) 

Dry 

Thick 

.62 

.38 

14(2} 

1187(.17) 

629(.09) 

2717(389) 

1544(221) 

1837(263) 

1048(150) 

1055(151) 

1055(151) 

Wet 

Thick 

.62 

.38 

14(2) 

978(.14) 

419(.06) 

2291(328) 

1055(151) 

1551(222) 

712(102) 

740(106) 

740(106) 

D-y 

Thin 

.66 

.34 

14(2) 

1187(.17) 

629(.09) 

2969(425) 

1697(243) 

2005(287) 

1153(165) 

1090(156) 

1090(156) 

Wet 

Thin 

.66 

.34 

14(2) 

1048(.15) 

489(.07) 

2689(385) 

1348(193) 

1816(260) 

908(130) 

957(137) 

957(137) 

Matrix  tangent  shear  modulus  corresponding  to  appropriate  matrix  shear  stress  from  analysis 
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Figure  3  Sandwich  Beam  Test  Specimen 

Copyright  ASTM.  Reprinted  with  permission. 


Core  material 


Test  face 
-L. 


X 


Opposite  face 


Notes. 

1 .  Load  pads  —  f  wide 

2 .  Reaction  pads  —  g  wide 


British 

Metric. 

system,  In. 

mm 

a 

4.0 

101  6 

b 

80 

203  2 

c 

1 .0 

25  40 

d 

200 

506  0 

e 

to 

25  4 

f 

1  0 

25  4 

0 

1  5 

38  t 
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Figure  4  Compression  Test  Method 
for  Rigid  Plastics 
ASTM  D695 
Copyright  ASTM. 
Reprinted  with  permission. 


Figure  5  ASTM  D695  Test  Specimen  Geometry 

Copyright  ASTM.  Reprinted  with  permission. 
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Figure  6  DTRC  Compression  Fixture 
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Figure  7  Failure  Modes  Corresponding  to  Highest 
Compression  Strengths  in  Unidirectional 
AS4/3501-6  Composites  Tested  Using 
the  IITRI  Compression  Fixture  [23] 


Figure  8  Kink  Band  Geometry 
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Figure  10  Orientation  Andes  for  Calculation 
Of  Moisture  Diffusion  In  A 
Composite  Material 


y 


TEMPERATURE  (°F) 


Figure  1 1  Curing  Schedule  for  AS4/3501-6  Panels 


TIME  (MIN) 


1)  Apply  28  inches  vacuum. 

2)  Heat  part  to  210°F  at  1-5°/min. 

3)  Hold  at  210°F  for  1  hour. 

4)  Apply  90  psi  pressure.  Vent  vacuum  when  pressure  reaches  25  psi. 

5)  Raise  temperature  to  355°F  at  1-5°/min. 

6)  Hold  temperature  at  355°F  for  2  hours. 

7)  Cool  to  150°F  at  rate  of  1-5°/min. 

8)  Vent  pressure  at  150°F. 

9)  Hold  at  355°F  for  4  hours. 

10)  Cool  to  75°F  at  rate  of  2 °/min.  and  vent  pressure. 
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Figure  12  Compression  Specimens 


Thick  Thin 

(All  dimensions  in  millimeters) 

Surfaces  A  and  B  ground  flat  and  parallel  to  0.03  mm 
Thickness  as  received  following  hand  layup  and  autoclave  cure 
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Figure  13  Method  Used  to  Determine  Coefficient 
of  Friction  Between  Specimen 
and  Fixture 

With  and  Without  Sandpaper 
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Figure  14  Elevated  Temperature  Setup 
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Figure  1 5  Schematic  of  DTRC  Thick-Section 
Compression  Fixture  After  [9] 


STRESS  (MPA) 


Figure  1 6  Compression  Stress-Strain  Curves 
for  Thick  Specimens  in  the 
Dry  Condition 


STRAIN  X  10E-6 
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STRESS 


Figure  17  Compression  Stress-Strain  Curves 
for  Thick  Specimens  in  the 
Wet  Condition 


-20000  -16000  -12000  -8000  -4000  0 

STRAIN  X  10E-6 
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STRESS 


STRESS (MPA) 


Figure  19  Compression  Stress-Strain  Curves 
for  Thin  Specimens  in  the 
Wet  Condition 


STRAIN  XI OE-6 
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Figure  20  Typical  Failure  Modes  for  Dry  and  Wet 

AS4/3501-6  (0  /90)  Thick  (6.35  mm)  Laminates 

L.  O 


9S 


Figure  21  Typical  Failure  Modes  for  Dry  and  Wet 

AS4/3501-6  (02/90)s  Thick  (6.35  mm)  Laminates 
Detail  View  of  Failure  Zone 


:9 


Figure  22 


Typical  Failure  Modes  for  Dry  and  Wet 
AS4/3501-6  (CL  /90V  Thin  (1.6  mm)  Laminates 

4  o 


DRY  WET 
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Figure  23  Typical  Failure  Modes  for  Dry  and  Wet 

AS4/3501-6  (02/90)sThin  (1.6  mm)  Laminates 
Detail  View  of  Failure  Zone 


DRY  WET 
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Figure  26  Boundary  and  Loading  Conditions 
for  3-D  Finite  Element  Analysis  of 
Compression  Specimens 


Surface  A 


Figure  27  Distribution  of  Moisture  Concentration 

Through  the  Thickness  for  Thick  (6.35  mm) 
AS4/3501  -6  (02/90)s  Laminates 


Figure  28  Distribution  of  Moisture  Concentration 

Through  the  Thickness  for  Thin  (1 .6  mm) 
AS4/3501-6  (02/90)s  Laminates 


Figure  29  Micromechanical  Finite  Element 
Model  of  Fiber  Plate  and 
Surrounding  Matirx  Layers 


Section  A  Section  A 


a)  Single  Fiber  Plate  Model 


Section  A 


Section  A 


b)  Double  Fiber  Plate  Model 
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SHEAR  STRESS  (MPA) 


Figure  30  3501-6  Resin  Shear  Stress-Strain  Curves 
Generated  with  Richard-Blacklock 
Shear  Model 
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UNITS  FROM  FIBER-MATRIX  INTERFACE 


Figure  31  Rotated  Shear  Stresses  at  the  Integration  Points 
on  Either  Side  of  the  Fiber-Matrix  Interface 
a/L  =  0.00375 


B  Gm=  1.60  GPA 
•  G;=  0.99  GPA 
■  G;=  0.37  GPA 


E  2 
=1 

CM 

<D 

r-  o 


FIBER 


SHEAR  STRESS  (MPA) 


Value  of  matrix  shear  modulus  used  in  finite  element  analysis 
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UNITS  FROM  FIBER-MATRIX  INTERFACE 


Figure  32  Rotated  Shear  Stresses  at  the  Integration  Points 
on  Either  Side  of  the  Fiber-Matrix  Interface 
a/L  =  0.009 


Value  of  matrix  shear  modulus  used  in  finite  element  analysis 
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